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Abstract: We update and extend to larger masses our previous analysis of the MSSM 
with minimal SOiq [MSOiqSM] soft SUSY breaking boundary conditions. We find a well- 
defined, narrow region of parameter space which provides the observed relic density of dark 
5^ \ matter, in a domain selected to fit precision electroweak data, including top, bottom and 

tau masses. The model is highly constrained which allows us to make several predictions. 
We find the light Higgs mass rrih < 121 =b 3 GeV and also upper bounds on the mass of the 
gluino m-g < 3.1 TeV and lightest neutralino rriy- < 450 GeV. As the CP odd Higgs mass 
niA increases, the region of parameter space consistent with WMAP data is forced to larger 
values of smaller values of m/^. Hence, we find an upper bound ruA < 1.3 TeV. 

This in turn leads to lower bounds on BR(i?s fi'^ jjL^) > 10^^ (assuming minimal flavor 
violation) and on the dark matter spin independent detection cross section a^^ > 10^^ pb. 
Finally, we extend our previous analysis to include WIMP signals in indirect detection and 
find prospects for WIMP detection generally much less promising than in direct WIMP 
searches. 



Keywords: Supersymmetric Effective Theories, Cosmology of Theories beyond the 



SM, Dark Matter. 



Contents 



|l]. Introduction ^ 

g. Minimal SOm SUSY Model MSOiqSM g 
|2.l| Framework ^ 



2.2 Phenomenological Analysis 



§. Results § 
|3.1| Experimental Tests ^ 
|3.2| Some representative points 19 



El Predictions and Summary 19 



1. Introduction 

The constrained minimal supersymmetric standard model [CMSSM] |Q] is a well defined 
model for soft SUSY breaking with five independent parameters given by rriQ, M1/2, Aq, 
tan (3 and sign{fj,). It has been used extensively for benchmark points for collider searches, 
as well as for astrophysical and dark matter analyses. The economy of parameters in 
this scheme makes it a useful tool for exploring SUSY phenomena. However the CMSSM 
misses regions of soft SUSY breaking parameter space which give qualitatively different 
predictions. In a previous paper Q we considered an alternate scheme, the minimal SOio 
supersymmetric model [MSOiqSM] 1^], which is well motivated and opens up a qualitatively 
new region of parameter space. In light of the recent WMAP analysis and recent improved 
limits from DZero and CDF on the branching ratio BR(i?s from CDMS on 

the dark matter direct detection cross section, we have decided to reanalyze the MSOiqSM, 
including updated data and extending the region of parameter space to larger values of 
mid > 3 TeV and niA > 500 GeV. 

The recent WMAP data [Q] provides an important constraint on the model. The dark 
matter candidate in this model is the lightest neutralino. However, since the scalar masses 
of the first two families are of order mig > 1.2 TeV, and the third generation sfermions 
(except for the stops) also tend to be heavy, the usually dominant annihilation channels, 
for the neutralino LSP to light fermions via t-channel sfermion exchange, are suppressed. 
On the other hand, the process XX ~^ ff via s-channel CP odd Higgs A exchange becomes 
important. This is due to the enhanced A coupling to down-type fermions, which is 
proportional to tan /?, and because, in contrast to heavy scalar exchange, the process is not 
p-wave suppressed. 
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We also compute the branching ratio for the process Bs —>■ due to A exchange ||5|, 

|6|. It is absolutely essential to include this latter constraint in our analysis, particularly 
in light of the published DZero bound BR{Bs /i+ /i") < 5.0 x 10"'^(4.1 x lO""^) at 95% 
(90%) CL and the new prehminary CDF bound < 2.0 x lO""^ at 95% CL [| and DZero 
bound < 3.7 x lO""^ at 95% CL §. 

The CDMS Collaboration has recently improved the upper limit on the spin inde- 
pendent dark matter WIMP elastic scattering cross section on a proton, cf^^ , down to 
2 X 10~^ pb (at low WIMP mass) [^]. A further improvement by an order of magnitude 
is foreseen within a year. We update our previous results for a^^ and show that a large 
fraction of the parameter space will be probed with currently running detectors, and will 
be completely explored with a new round of "one-tonne" detectors which plan to reach 
down to a^^ > 10~^° pb. 

In addition, we now evaluate prospects for indirect detection of WIMPs in muon flux 
from WIMP annihilation in the Sun, a gamma ray flux from the Galactic center, as well 
as antiproton and positron fluxes from the Galactic halo. Here detection prospects are 
somewhat varied and, especially in the case of the gamma ray flux, strongly depend on an 
adopted model of the Galactic halo. 

The paper is organized as follows. In section |2| we review the MSOiqSM, describe its 
virtues and outline the analysis. Then in section ^ for different values of mig and rriA we 
compute BR(i?s — > fj.'^ M~)j the cosmological dark matter density fi^/i^, etc. The major 
results of the paper are found in section 3.1. We discuss the upper bound on the CP 
odd Higgs mass, mA_, and the resulting lower bound on BR(i?s /i+ fi^). CDF and 
DZero can potentially probe the entire allowed range. We also consider the predictions 
for underground dark matter searches. The entire expected range can be fully explored at 
the next level of underground dark matter searches. We also consider indirect dark matter 
searches. Finally in Table 1 we give Higgs and SUSY spectra relevant for collider searches 
for four representative points in SUSY parameter space. We summarize our results in 
section 



2. Minimal SOm SUSY Model - MSOiqSM 
2.1 Framework 

Since there are several distinct theories in the literature called the minimal SO(IO) SUSY 
model, let us briefly describe here the properties of our minimal SOiq SUSY model 
[MSOiqSM] |^]. Quarks and leptons of one family reside in the 16 dimensional repre- 
sentation, while the two Higgs doublets of the MSSM reside in one 10 dimensional repre- 
sentation. For the third generation we assume the minimal Yukawa coupling term given by 
A 16 10 16. On the other hand, for the first two generations and for their mixing with the 
third, we assume a hierarchical mass matrix structure due to effective higher dimensional 
operators. Hence the third generation Yukawa couplings satisfy At = = = Xu^ = 
Soft SUSY breaking parameters are also consistent with SOiq with (1) a universal 
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gaugino mass Mx/2) (2) a universal squark and slepton mass miQ,^ (3) a universal scalar 
Higgs mass mio, and (4) a universal A parameter Aq. In addition we have the supersym- 
metric (soft SUSY breaking) Higgs mass parameters fj, {B^). Bfi may, as in the CMSSM, 
be exchanged for tan /?. Note, not all of these parameters are independent. Indeed, in order 
to fit the low energy electroweak data, including the third generation fermion masses, it 
has been shown that Aq, mio, miQ must satisfy the constraints Q 

^0 ~ -2 mie; mio ~ niiQ (2.1) 

mi6 > 1.2 TeV; //, M1/2 <C rriie (2.2) 

with 

tan/3«50. (2.3) 



This result has been confirmed by several independent analyses [11, |T2|.^ Although the 
conditions (Eqns. |2.2D are not obvious, it is however easy to see that (Eqn. (|2.3| )) 
is simply a consequence of third generation Yukawa unification, since mt{mt) / m}j{mt) ~ 
tan j3. 

One loop threshold corrections at the GUT scale lead to two significant parameters we 
treat as free parameters, although they are calculable in any GUT. The first is a correction 
to gauge coupling unification given by 

eg = [a^iMc) - ac] /ac (2.4) 

where the GUT scale Mq is defined as the scale where ai(MG') = a2{MG) = do- The 
second is a Higgs splitting mass parameter defined by 

Ami = {ml, - mlj/2mlQ. (2.5) 

In order to fit the low energy data we find £3 ~ —4% and Amfj ~ 13% Q. The largest 
corrections to €3 come from the Higgs and SOiq breaking sectors, while the correction to 
Ami predominantly due to the right-handed r neutrino. Note, for Mp^ ~ 5.8 x 10^^ 
GeV, the necessary Higgs splitting is reduced to 7% |14|. 



Finally, as a bonus, these same values of soft SUSY breaking parameters, with mie 'S> 
TeV, result in two very interesting consequences. Firstly, it "naturally" produces an in- 
verted scalar mass hierarchy [ISMH] [^]. With an ISMH squarks and sleptons of the first 
two generations obtain mass of order miQ at Mz- The stop, sbottom, and stau, on the other 
hand, have mass less than (or of order) a TeV. An ISMH has two virtues. (1) It preserves 
"naturalness" (for values of mig which are not too large), since only the third generation 
squarks and sleptons couple strongly to the Higgs. (2) It ameliorates the SUSY CP and 
fiavor problems, since these constraints on CP violating angles or fiavor violating squark 
and slepton masses are strongest for the first two generations, yet they are suppressed as 



^SOio does not require all sfermions to have the same mass. This however may be enforced by non- 
abelian family symmetries or possibly by the SUSY breaking mechanism. 

^Note, different regions of parameter space consistent with Yukawa unification have also been discussed 
in[0|l||l|l. 
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1/mfg. For mi6 > a few TeV, these constraints are weakened |T^, 17, 18 1. Secondly, Super- 
Kamiokande bounds on t{p K^v) > 2.3 x 10^^ yrs |1£] constrain the contribution of 
dimension 5 baryon and lepton number violating operators. These are however minimized 
with //, Ml/2 < "^16 

2.2 Phenomenological Analysis 

We use a top-down approach with a global analysis [21|. The input parameters are 
defined by boundary conditions at the GUT scale. The 11 input parameters at Mq are given 
by: three gauge parameters Mq, adMo), £3; the Yukawa coupling A, and 7 parameters, 
including pi, and the 6 soft SUSY breaking parameters M1/2, ^0, tan/3 (replacing B^), 
m^g, m^Q, Am'jj. These are fit in a global analysis defined in terms of physical low 
energy observables. Note we keep three parameters (mig, fi, M1/2) fixed; while minimizing 
with the remaining 8 parameters. Below we plot contours as a function of fi, Mi^2 
for different values of mig. We use two (one) loop renormalization group [RG] running 
for dimensionless (dimensionful) parameters from Mq to Mz- ^ We require electroweak 
symmetry breaking using an improved Higgs potential, including mf and mf corrections 
in an effective 2-Higgs doublet model below MsusY = \J\{'''^\ + lH' 

The function includes 9 observables; 6 precision electroweak data a^M, G^, as{Mz), 
Mz, Mw, PNEW and the 3 fermion masses Mtop, mh{mb), Mr- In our analysis we fit 
the central values pi: Mz = 91.188 GeV, Mw = 80.419 GeV, x 105 = 1.1664 



GeV ^, = 137.04, M^- = 1.7770 GeV with 0.1% numerical uncertainties; and the 

following with the experimental uncertainty in parentheses: as{Mz) = 0.1172 (0.0020), 
Pnew X 103 = -0.200 (1.1) H], Mt = 178.0 (4.3) GeV, mb{mb) = 4.20 (0.20) GeV. ^ We 
include the complete one loop threshold corrections at Mz to all observables. In addition 
we use one loop QED and three loop QCD RG running below Mz- 

The output of this analysis is a set of weak scale squark, slepton, gaugino and Higgs 
masses. With regards to the calculated Higgs and sparticle masses, the neutral Higgs 
masses h, H, A are pole masses calculated with the leading top, bottom, stop, sbottom 
loop contributions; while all other sparticle masses are running masses. This output is then 
used to compute the cosmological dark matter density of the lightest neutralino Q.^h?, which 
is the LSP, the branching ratio BR(5s — > p^ P~), and other observables, as described in 
more detail below. Note, it is important to emphasize that neither r^^^/i^ nor the branching 
ratio BR(i?s — > p'^ p~) are included in the x'^ analysis. They are predictions of the model 
obtained in the regions selected by consistency with gauge coupling and third generation 
Yukawa unification and the low energy observables. 

Using x^ penalties^ we apply two additional constraints: 



^Note, we have checked that switching to 2 loop RGEs for dimensionful parameters can be compensated 
for by small changes in the GUT scale parameters, without significant changes in the low energy results. 

■^Note we take a conservative error for mtimt) in view of recent claims to much smaller error bars [ p6[ . 

^In order to constrain the values of some physical observables in our analysis, such as m^^ or tua, we 
add a significant contribution to the function for values of these observables outside the desired range. 
We refer to this additional contribution as a penalty. Minimization of x^ with Minuit, then pushes the 
fits to the desired range. Of course the x^ penalties then vanish. 



- 4 - 



• > 300 GeV 

• niA fixed. 

The first is cliosen to be consistent with BR(i? Xg'y) |^]. Note, although we do 
calculate BR(i3 ^sT)) we do not use it as a constraint in the analysis. This is for two 
reasons — 1) this decay mode depends on 3-2 generation mixing which is model dependent 
and 2) it is not difficult to fit BR(i? ^sl) for large enough values of m^^. Hence, in 
order to be generally consistent with the measured value of BR(i? — > Xg^)^ we impose 
TTiii — GeV. With regards to the second constraint, since ^i^/i^ and BR(i?s — > /i^ 
are both sensitive to the value of niA, we fix its value and present our results for different 
values of tua-^ 

Finally, after performing the analysis, we impose the following constraints on the 
parameter space: 

• lower bound on the lightest chargino mass ?tT';^+ > 104 GeV. 

• lower bound on the light Higgs mass m/j > 111 GeV. Note, because of the theoretical 
uncertainty in the calculation of m/j (~ 3 GeV), we conservatively impose rrih > 
111 GeV, instead of the LEP bound for SM Higgs ruh > 114.4 GeV. 

• branching ratio BR(i?s /i+ /J^)- The current best published limit comes from 
DZero BR(5, ^ /i+ /i") < 5.0 x 10~'^(4.1 x IQ-"^) at 95% (90%) CL §. Recently 
CDF has announced a new preliminary bound < 2.0 x 10^^ at 95% CL |Q, while 
DZero has come up with a new preliminary bound < 3.7 x lO"'^ at 95% CL g. We'll 
discuss the impact of the bounds below. 

• the relic abundance of the lightest neutralino in the range 0.094 < il^/i^ < 0.129 (2 a) 
which we take as a preferred range. We will exclude points for which J^^^^ ^ 0.129 
but allow for the possibility that Q^h'^ < 0.094 (subdominant component of cold 
DM). 

The quantity BR(i?s ^'^ is computed assuming the CKM mixings among 
squarks (minimal flavor violation). For this purpose we use leading log expressions derived 
in the first paper of Ref. |^]. The SM prediction is around 3 x 10~^ while SUSY contribu- 
tions to the process are dominated by pseudoscalar exchange and scale as tan^ f3/m\, and 
can be large. Note, the branching ratio also depends on the (model dependent) 3-2 genera- 
tion squark mixings. However, the mixings tend to increase BR(i?s fi"^ /i~) which make 
our lower bound on rriA stronger. On the other hand, beyond leading log corrections can 



relax our bounds on ruA by up to 0(20%) because of a "focusing effect" |17, 18], especially 
if the squark mixings are non-minimal. 

We compute the relic abundance O^/i^ of the lightest neutralino using exact expressions 
for neutralino pair annihilation into all allowed final-state channels, which are valid both 



®The calculation of BR(i3 ^s7) and BR(i3i, /i+ fj,~) requires a model for fermion mass matrices. 
In the absence of such a model we use the observed CKM matrix elements to calculate these flavor violating 
branching ratios. 
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Figure 1: Contours of constant for mig = 3 TeV and vtla = 500 GeV (left panel) and 700 GeV 
(right panel). The red (light shaded) regions are excluded by < 104 GeV (below and to the 
left of the solid black curve), m/j < 111 GeV (on the right) and by fi^/i^ > 0.129. To the right of 
the thick broken black line one has trih < 114.4 GeV. The green (darkest shaded) band corresponds 
to the preferred 2 a range 0.094 < ^l^h^ < 0.129, while the white regions below it correspond to 
ri^/i^ < 0.094. The region excluded by the DZero experimental bound on BK{Bs ^ /i+ /i~) < 
5.0 X 10"^ at 95% CL is marked in dark blue (dark shaded), while the region affected by the new 
preliminary CDF bound < 2.0 x 10^^ at 95% is marked in light blue (light shaded). Contours of 
constant BR(i?s /i~) are given by the blue dashed lines. 



near and further away from resonances and thresholds [^] . We further treat the neutralino 
coannihilation with the lightest chargino and next -to-lightest neutralino |28] and with the 
lighter stau with similar precision. We include all coannihilation channels, including 
the previously neglected Q coannihilation with light stops, although this channel only 
affects points in the parameter space which are excluded by other constraints. We solve 
the Boltzmann equation numerically as in [^] and compute ri^/i^ with an error of a few 
per cent, which is comparable with today's accuracy on the observational side. 



3. Results 



In Figs. ^, |2| and ^ we present our results for different values of mig and tua in the fi, M1/2 
plane. For example, in Fig. |^ we present, for mie = 3 TeV and niA = 500 GeV (left 
panel) and niA = 700 GeV (right panel) the (magenta) solid lines of constant x^- The red 
(lightest shaded) regions are excluded by collider limits and by J^^^^ > 0.129. Specifically 
the red regions at small fj, and/or small M1/2 (bounded by the solid black line) is excluded 
by bounds on the chargino mass > 104 GeV and the red region on the right side 
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Figure 2: Same as Fig. | for mig = 5 TeV and niA = 700 GeV (left window) and 1250 GeV 
(right window). 



(bounded by the dashed black line) is excluded by (our conservative application of) the 
LEP Higgs bound ruh > 111 GeV. The cosmologically preferred dark matter region green 
(darkest shaded) satisfy 0.094 < ^^h^ < 0.129. We find significant regions of parameter 
space which give ^ 3, satisfy the preferred ^x^^ range as above, and satisfy all other 
phenomenological constraints. In Fig. ^, we have mig = 5 TeV and = 700 and 
1250 GeV, respectively, and in Fig. ^ we have mig = 3 and 5 TeV, respectively, and 
mA = l TeV. 

In Fig. H, marked in dark blue (dark shaded), is the region excluded by the DZero 
experimental bound on BR(i?s fi+ /i") < 5.0 x 10""^ at 95% CL 0] and, in Figs. | 
and ^, light blue (light shaded) is the region affected by the new preliminary CDF bound 
< 2.0 X 10"'' at 95% CL |^]. In addition, we have included contours of constant BR(ii?s 
/i+ H^) with dashed blue lines. The branching ratio BR(i?5 — > fi^) is sensitive to the 
value of the CP odd Higgs mass rriA scaling as m'^^. For niA = 500 GeV (Fig. || (left)) 
the branching ratio BR(i?s fi^ fi^) is below the published DZero bound, for acceptable 
values of ^x^"^ ^ 3, but is almost excluded by the preliminary CDF bound. While 

for rriA = 1-25 TeV (see Fig. § (right window)) we have BR(i?s fi^ fi^) > 10^^. 

The cosmological relic abundance of the neutralino f^^^^ primarily determined by the 
direct s-channel pair-annihilation into SM fermion pairs through the CP odd Higgs. Since 
all the sfermions are very heavy, their contribution to reducing the neutralino number 
density is strongly suppressed. In contrast, because of the coupling Abb oc tan (3 (and 
similarly for the r's), the ^-resonance is effective and broad. Near nix ~ ™'A/2 it reduces 
Qxh'^ down to allowed but uninterestingly small values <C 0.1. As one moves away from 
the resonance, ^i^/i-^ grows, reaches the preferred range 0.094 < fi^^/i'^ < 0.129, before 
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Ml/, (GeV) Mj/, (GeV) 

Figure 3: Same as Fig. for mig — 3 TeV (left window) 5 TeV (riglit window) and rriA = 1 TeV. 



becoming too large ^x^^ 0.129. ^ (A similar, but much more narrow resonance due to 
hP is also present at M1/2 ~ 150 GeV and small /i.) When > rrit (M1/2 ^ 420 GeV) 
and the stops are not too heavy, the LSP pairs annihilate to t t-pairs. In the region of large 
Ml/2, often where mh is already too low, two additional channels become effective. First, 
in this region the neutralino becomes almost mass degenerate with the lighter stau which 
leads to reducing J^^^/i^ through coannihilation. Second, if tua is not too large, neutralino 
pair-annihilation into Higgs boson pairs AA and HH opens up. Finally, at <C M1/2, the 
relic abundance is strongly reduced due to the increasing higgsino component of the LSP. 

3.1 Experimental Tests 

BR{Bs — >■ /x^/x"). As we see in Figs. |l|-^, the preliminary Tevatron limit on BR{Bs — > 
Iji^ lJi~) puts a significant lower bound on > 500 GeV. The dependence on niA is plotted 
in Fig. ^ with points in the preferred region with all the collider constraints satisfied, \^ <2> 
and 0.094 < Q^/i^ < 0.129, for the two representative cases mig = 3 and 5 TeV. Note that, 
for a fixed value of m^, a vertical spread of points is caused by the changing shape of the 
cosmologically favored band in the /i, Mij2 plane and some variation in tan j3 in Figs. ^[-|3|. 

We now argue that there is an upper bound on ruA, the CP odd Higgs mass. Hence 
there is a lower bound on the branching ratio BR{Bs — > We first show that 

increasing the value of mig permits a larger range for the parameters fi, Mi/2- Comparing 
Figs. |l|-|3| we see that, as mig increases, the region with < 3 rapidly grows. Note, 
the dominant pull in is due to the bottom quark mass. In order to fit the data, the 
total SUSY corrections to mi,(mi,) must be of order —(2 — 4)% In addition there 

'^Note that at one loop we have Mi(Mz) = M1/2 * ai{Mz)/aG so Mi(Mz) ^ 0.4Mi/2. For bino-Uke 
neutraUno (which is true for larger jj,), we thus have ~ 0.4A/i/2- Hence for s-channel annihilation we 
have niA ~ 2m^ « 0.8 M1/2 or M1/2 ~ (5/4) niA for the position of the "peak suppression." 
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Figure 4: The branching ratio BR{Bs /^^) as a function of in the /i, M1/2 region of 

parameter space satisfying all the collider constraints, 0.094 < Q,^h? < 0.129 and < 3 for fixed 
mi6 = 3 and 5 TeV. 



are three dominant contributions to these SUSY corrections: a gluino loop contribution 
oc Q!3 /i Mg tan/3/m? , a chargino loop contribution oc fi At tan/?/m?^, and a term 
oc log Mg^gY■ Larger values of uiiq permit a larger range for the ratio m^^/m-j^. Thus 
larger values of mig allows more freedom in parameter space for fitting the data at both 
smaller or larger values of Mx/2- Note, also when mig increases (with M1/2 fixed) the 
parameter At becomes more negative, since Aq ~ —2 niiQ and At ~ —3 M1/2 + ^ ^0 where 
e < 1. 



In addition, as increases the light Higgs mass decreases. The reason for this can 
easily be seen. Consider the approximate equation for the light Higgs mass (Eqn. (2.32) 
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from Ref. 



, 3 mt{mt)^ 



47r2 



lX, + t + j^ (i^^^ - 32^a3) {Xtt + t^)] (3.1 



where Xt = .1^* (l - ttjTt ), t = \og{Ml^ gy / ^t{mtf) and MsusY 



^3uSY \ 12AfJ 

(m? + m? )/2. In our analysis, the bottom quark gets a positive supersymmetric correc- 
tion proportional to and a negative contribution proportional to ^At. Thus in order 

to fit the bottom quark mass as increases, —At also increases. As a consequence, in 

the relevant region of parameter space, Xt decreases. In addition, the fit value of mt{mt) 
also decreases. Hence the value of rrih decreases as Mi 12 increases. On the other hand, 
the value of the light Higgs mass is fairly insensitive to rriig.^ In the acceptable regions of 
parameter space we find 114.4 < ruh < 121 GeV. ^'^ 

Now let us put it all together. Comparing various windows in Figs. |l|-^, we see that 
increasing tua has two effects. It suppresses the branching fraction BR(i?s — > /^~)- 
At the same time it moves the s-channel neutralino annihilation channel to larger values 
of Ml providing larger regions with 0.094 < VLy^h'^ < 0.129. However, mcreasmg 
above 1 TeV or so moves the regions of preferred ^^^^^ ^° right, in potential 

conflict with a lower bound on ruh- In fact, as can be seen from Fig. || (1 ower right), 
there is a rough upper bound on of about 1.3 TeV {and a corresponding lower bound on 
BR(5s —>■ fi^ > O{10~^)) , above which there are no longer any solutions consistent 
with the observed fi^^/i^ and the lower bound on the Higgs mass. 

It has been shown that, with an integrated luminosity of 15 fb~^, experiments at the 



Tevatron are sensitive to BR(i?s /U"*" > 1.2 x 10^^ |31]. Hence this offers the 
prospect of a full exploration of the model in this observable. In addition, the LHC is 
expected to be sensitive down to (3.5 it 1.0) x 10~^ |32|. We discuss BR{Bs — > /U^/i^) 
further below. 

WIMP Direct Detection Search. We begin by plotting in Fig. |5| a^^ , the spin in- 
dependent neutralino dark matter cross section relevant for direct dark matter searches, 
versus tua- We show the cases mig = 3 TeV, tua = 500, 700, 1000 GeV (dark blue 
crosses) and mig = 5 TeV, tua = 500, 700, 1000, 1250 GeV (magenta circles). The 
points lie in the preferred region with all the collider constraints satisfied, < ?> and 
0.094 < < 0.129. Like in Fig. |, for a fixed value of niA, a vertical spread of points is 

caused by the changing shape of the cosmologically favored band in the M112 plane and 



*This approximation gives a value for the light Higgs mass which is larger than the value obtained with 
a more exact calculation (see Fig. 2 in Ref. for values of \At \ above 2 TeV. The discrepancy increases 
with increasing |j4t| > 2 TeV, overestimating the Higgs mass by as much as 8 GeV for \At \ = 3.6 TeV. 

^Although for fixed M-1/2 the light Higgs mass does increase by a small amount as mie increases (see 
Fig-i). 

^"Note, there is at least a ±3 GeV theoretical uncertainty in the Higgs mass. 

^^The upper bound on tua will increase slightly as mig increases, but then mie cannot increase signifi- 
cantly without creating a fine-tuning problem. 
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Figure 5: Spin independent neutralino dark matter cross section versus niA for miQ — 3 TeV 
(dark blue crosses) and 5 TeV (magenta circles). The points obey all collider constraints, < 3, 
0.094 < n^h'^ < 0.129 and the new preliminary CDF bound BR{Bs -> fi+ fJ.^) < 2.0 x lO"'^. 



also by some variation in tan /? in Figs. ||-|3|. The dominant contribution to both tree-level 
and one-loop diagrams typically comes from a i-channel exchange of the heavier scalar 
Higgs boson H'^. Thus one would normally expect 0"^^ ~ "^a^ since ~ ttih- However, 
the coupling xxH^ is sensitive to the bino/higgsino composition of the lightest neutralino 
and vanishes in the limit of the neutralino becoming one of these two states. The bino 
fraction increases along the cosmologically favored band from as low as some 70% at low 
H to almost 100% at largest allowed jjL. This causes a^^ to decrease with ji even for a fixed 
ruA- 

Since both BR{Bs fJ-^f^ ) and are to a large degree determined by ruA, but 
also by the LSP composition, one finds interesting correlations which are shown in Fig. |6|. 
The small light points obey all collider constraints and < 3, while the larger dark 
points satisfy in addition 0.094 < J^^/i^ < 0.129. For fixed tua and M1/2, cf^^ decreases 
with increasing /i due to the increasing bino content of the neutralino. On the other hand 
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Figure 6: The branching ratio BR{Bs ^ /j,"*" versus spin independent neutralino dark matter 
cross section a^^ for niie = 3 TeV (left window) and 5 TeV (right window) and for several choices 
of TJiA- The small light points obey all collider constraints and < 3, while the larger dark points 
satisfy in addition 0.094 < fl^h'^ < 0.129. 



BR(i?s fi'^ n^) increases as fi increases, since the flavor violating vertex correction 
is proportional to fi (see below). Hence, as seen from Fig. |6|, in the allowed parameter 
space of the MSOioSM, the values for (7^^ and BR(i?s — > fj,'^ fi^) are (as one moves 
along the cosmologically favored band, for a fixed m^) inversely correlated. This behavior 
sharply contrasts with the direct correlation between the two quantities in the CMSSM 
and its extensions presented in |3S], where, with the usual parametrization, the masses of 
the Higgs bosons grow along the cosmologically favored band. Thus, an improved limit on 
BR{Bs — > does not need to necessarily imply less promising prospects for measuring 

(Tp^ , and vice versa. In addition, in Fig. |6| we see that, comparing the case of mig = 3 and 
5 TeV (again, for fixed niA) in the latter case the range of a^^ for the (dark) points in 
the cosmologically preferred region extends to larger values of BR(i3s l^~)- This is 

because the flavor violating vertex of the CP odd Higgs scales roughly like //A^/m^g and 
for larger values of mig the acceptable regions of parameter space extend to even larger 
values of ji. 

In Fig. 1^ we plot a^^ versus the neutralino mass. We show the same cases of mig and 
niA as in Figs. Q-^. Big points obey all collider constraints, < 3, 0.094 < ^^h^ < 0.129 
and the new preliminary CDF bound BR(i?s fj.^ fi^) < 2.0 x 10^^ while for the small, 
light points the last constraint is relaxed to BR(i?s — > fi^ jj,^) < 5.0 x 10"^. In the (light 
blue) background we plot predictions obtained in the general MSSM by scanning large 
ranges of parameters, including tan/3 from 10 up to 65 |34] but do not apply the bound 
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Figure 7: Spin independent dark matter cross section tip ^ vs. the neutralino mass for several 
choices of mie and niA- Big points obey all collider constraints, < 3, 0.094 < il^h^ < 0.129 and 
the new preliminary CDF bound BR(i?s ji^ < 2.0 x 10^^ while for the small, light points 
the last constraint is relaxed to BR(i?s ^ /i+ /i^) < 5.0 x 10^^. The hghtly shaded background 
represents predictions in the general MSSM for which the bound from BR(_Bs /i+ was not 
applied. 



from BR(i?s /^^ ^ )• Note, most regions of parameter space should be observable 
as dark matter search experiments are expected to probe down to values of a^^ > 10^^ 
pb with current /upgraded detectors within the next year or so {e.g., by CDMS). Future 
one-tonne detectors are planned to reach cjp^ > 10^^'^ pb, thus probing the entire favored 
parameter space of the MSOiqSM. 

WIMP Indirect Detection Search. Apart from direct detection, dark matter WIMP 
signals can be detected in a number of indirect detection search channels. Here we com- 
pute detection rates for a high-energy neutrino flux from the Sun, for gamma rays from 
the Galactic center, and for antiprotons and positrons from the Galactic halo, by applying 
DARKSUSY [^] to some popular Galactic halo models. As we will see, prospects for 
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Figure 8: Muon flux from the Sun vs. the neutrahno mass for tua = 1 TeV and mig = 3 TeV 
(left window) and mig = 5 TeV (right window) assuming the energy threshold Ejj^ = 25 GeV. 
In light blue (shaded) areas the collider constraints and < 3 have been applied but not the 
cosmological constraint 0.094 < fi^/i^ < 0.129 nor the bound on BR(i3s /i^ )• Points along 
the dark blue line additionally satisfy the cosmologically preferred range 0.094 < fl^h'^ < 0.129. 
The current limit from SuperKamiokande is marked as a red solid line while future reach of IceCube 
is marked as a black dashed line. 



WIMP detection in these modes are generally less promising than in direct detection ex- 
periments, although a sizeable number of cosmologically preferred cases can still be probed. 
On top of this, astrophysical uncertainties are often sizeable, especially for the latter two 
modes since propagation of antiprotons and positrons in the Milky Way is poorly under- 
stood. 

Instead of undertaking a full presentation of all the cases discussed above, we will select 
a representative value of niA = 1 TeV and will present our results assuming all collider 
constraints (except for the BR(i?s t-i-^ l^ ) bound which, for the chosen value of ruA, 
is irrelevant), < 3 and Q^^h? < 0.129, and then show the (typically strong) impact of 
imposing the cosmological constraint 0.094 < O.^h'^ < 0.129. Whenever the WIMP relic 
density is too small (Jl^/i^ < 0.094), we will not rescale the rates. 

Halo WIMPs pass through the Sun and other celestial bodies and occasionally become 
gravitationally trapped inside their cores. Once enough of them have accumulated, they 
start pair-annihilating into pairs of SM particles {qq, WW, ZZ, etc), which subsequently 
decay via two- and three-body processes, like b clvi, W vv, etc. Among the decay 
products only neutrinos can escape out of the Sun's core. A muon neutrino from the Sun, 
when passing through the Earth, may produce a muon which will generate a Cerenkov 
shower in a nearby under-ice or under-water detector. 
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Figure 9: Gamma ray fluxes from the Galactic center vs. the neutralino mass for = 1 TeV 
and for mig = 3 TeV (left window) and rriig = 5 TeV (right window), assuming the energy threshold 
= 1 GeV, as in GLAST. Results for three different halo models are shown, as described in the 
text. In light-colored (shaded) areas collider constraints and < 3 have been applied but not the 
cosmological constraint 0.094 < fi^/i^ < 0.129 nor the bound on BR{Bs fi~). Dark lines 

corresponding to each halo model mark points which additionally satisfy the cosmological constraint 
0.094 < n^h^ < 0.129. The current limit from EGRET is marked as a red solid curve while an 
expected reach of GLAST is marked as a black dashed line. 



The calculation of the WIMP capture rate and subsequent annihilation, as well as prop- 
agation of decay products is rather well understood and fairly insensitive to the choice 
of a Galactic halo model. (We assume the local dark matter density po = 0.3 GeV/cm^ and 
a Maxwellian velocity distribution with the peak velocity vq = 270 km/s.) The capture 
rate in the Sun, typically determined by spin dependent interactions involving Z-boson ex- 
change, is basically independent of m^. On the other hand, the shape of the cosmologically 
favored band in the fi, plane does depend on niA- 

In Fig. P we plot the high-energy muon flux from the Sun for uia = 1 TeV, and mig = 
3 TeV (left window) and mig = 5 TeV (right window). The IceCube energy threshold £'*^ = 
25 GeV has been applied. In light blue (shaded) areas the collider constraints and < 3 
have been applied but not the cosmological constraint 0.094 < Q^/i^ < 0.129 nor explicitly 
the bound on BR(i?s fi'^ fi^). Dark lines correspond to points which additionally satisfy 
the cosmologically preferred range 0.094 < ^l^h'^ < 0.129. We can see that the favored 
configurations fall below the current best limit from Super-Kamiokande [|l9| ] and are on a 
bordeline for being probed at IceCube [^] . 

Predictions for the fluxes of gamma rays from the Galactic center and antiprotons and 
(to a lesser extent) positrons from the Galactic halo depend on the assumed halo model. 
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Several profiles of dark matter distribution have been discussed in the literature, many of 
which can be parametrized by 



p{r) = po [1 + {ro/an^ , (3.2) 

[l + (r/a)"] — 

where p{r) is the radial dependence of the halo WIMP density, po = 0.3 GeV/cm^ is the 
local dark matter density, tq is our distance to the Galactic center and a is a distance scale. 
Here we consider three distinct and popular choices: 



a sperically symmetric modified isothermal profile [g3> for which (a,/?, 7) = (2,2,0), 
ro = 8.5 kpc and a = 3.5 kpc; 



the Navarro, Frenk and White (NFW) profile @, for which (a,/3,7) = (1,3,1) 
ro = 8.0 kpc and a = 20 kpc; 



• the Moore, et al, profile ||39|], for which (a,/3,7) = (1.5,3,1.5), ro = 8.0 kpc and 
a = 28 kpc. 

Towards the Galactic center the halo density is expected to be larger than in our 
local neighborhood but otherwise constant in the case of the modified isothermal model, 
or divergent as p(r) ~ 1/r (NFW) and p(r) ~ l/r^-^ (Moore, et al.). It is therefore natural 
to expect enhanced neutralino WIMP annihilations in the core of the Galactic center. 
Among annihilation products, high-energy gamma rays are unique in that they point back 
directly to the Galactic center. The flux is proportional to the DM number density squared 
integrated along the line of sight. Monochromatic photons with energy E.y ~ rriy. would 
provide a spectacular signal but the rates are low compared to the diffuse gamma radiation 
from cascade decays of WIMP annihilation products. 

In Fig. ^ we plot the diffuse high-energy gamma flux from the Galactic center, in- 
tegrated over the cone of 0.001 sr. We apply the continuous gamma energy threshold 
^th ^ I QgY^ planned for GLAST |o|. We also denote an upper hmit from EGRET |l| 



and the expected reach of GLAST after three years of taking data. The rates are indepen- 
dent of the choice of tjia- Depending on the choice of the halo model, GLAST may have a 
good chance of detecting a WIMP signal. 

For comparison, in Fig. ^ we plot the same quantity but apply the larger energy 
threshold of 60 GeV, typical of HESS. For the case mig = 5 TeV the cosmologically favored 
values of are higher than for m-ig = 3 TeV, and accordingly the resulting diffuse gamma 
radiation will be more accessible to HESS [42| than to GLAST. Thus HESS has a good 



chance of detecting a signal in the region of the parameter space which is to some extent 
complementary to the reach of GLAST. 

Another way of looking for WIMPs is to look for an excess of antiparticles in cosmic rays 



originating from WIMP annihilation in the Galactic halo, despite large uncertainties |35 



43]. A flux of antiprotons from astrophysical sources is expected to fall off dramatically at 



low energies, unlike that from WIMP annihilation in the halo. For a detailed discussion of 



antiproton fluxes from WIMP annihilation and from background sources see Ref. [43] 
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In Fig. 11 we plot a differential antiproton flux /dEdO, vs. m^. We evaluate the 
quantity assuming the peak of the flux at Ep- ~ 1.76 GeV to coincide with a peak in the 
kinetic energy distribution of the BESS experiment [44|. The plots do not depend on tjia 
and the halo model dependence is fairly weak, except for the Moore, et al, model for which 
the rate is somewhat higher. As we can see, the cosmologically favored points fall below 
the BESS sensitivity. 

Finally, high energy positrons may be a signature of WIMP annihilation in the halo. 
Again, like in the case of antiprotons, background from astrophysical processes (mostly 
spallation of cosmic rays off the interstellar medium) and other astrophysical uncertainties 
are rather large. Production of hard positrons in e~^e~ pairs from direct WIMP annihilation 
is tiny due to helicity suppression. On the other hand, annihilation into WW and ZZ pairs, 
followed by their decays into e~^e~ pairs will result in a positron flux peaked at an energy 
of ~ This feature gives some hope of distinguishing the signal (S) due to WIMP 

annihilation from the background (B) which is typically two to three orders of magnitude 
larger. 

In Fig. [l^ we plot the signal-to-background ratio S/B. The positron flux is evaluated 
at an optimal energy of and is compared with background flux, as parametrized by 

Feng, et al. |45|. Halo model dependence is negligible since energetic positrons produced 
in the solar neighborhood would be detectable. A sensitivity of positron search experi- 
ments PAMELA [46| and AMS-II [47| is indicated where S/B rates as low as 0.01 may be 
detectable. As can be seen from the figure, this is unlikely to be sufficient to probe our 
predicted rates for cosmologically favored regions. 
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Figure 10: Same as Fig. || but assuming the energy threshold E]^^ — 60 GeV, as in HESS whose 
limit is marked as a red solid line. 
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Figure 11: Same as Fig. g but for antiproton fluxes from the Galactic halo. A sensitivity of BESS 
is indicated with a dashed hne. 
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Figure 12: Same as Fig. ^ but for positron flux from WIMP annihilation in the Galactic halo 
presented as signal-to-background ratio B/S, as described in the text. Rates as low as 0.01 may 
be detectable at PAMELA and AMS-II. 
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3.2 Some representative points 



In Table 1 we present the input parameters and resulting Higgs and SUSY spectra for 
four representative points in SUSY parameter space which are consistent with all available 
data. We also present several other relevant experimental observables, including the SUSY 
contribution to the anomalous magnetic moment of the muon, a^^^ , and the branching 
ratio for the process B — > X^^, for these points. Note that one stop and/or one stau, and 
the charginos and neutralinos are relatively light. Also, as might be expected with heavy 
first and second generation sleptons, the value of af^^^ is small (of order 1-3 xlO"^*^). 
Finally, although BR(i? ^sT) was not included in x^, and the predicted values are 
typically outside the experimentally allowed range, it is not difficult to find consistent 
solutions to all the data, including this observable. 



4. Predictions and Summary 

This paper extends the analysis of a previous paper |§] to larger values of mig, from 3 to 
5 TeV, and larger CP odd Higgs mass ttla > 500 GeV. We find a well-defined, narrow 
region of parameter space which provides the observed relic density of dark matter, as 
well as a good fit to precision electroweak data, including top, bottom and tau masses, 
and acceptable bounds on the branching fraction of Bg ^ /i^ /i^. We present predictions 
for Higgs and SUSY spectra (Table), the dark matter detection cross section a^^ and the 
branching ratio BR(i?5 — > ^i^ in this region of parameter space. 

The MSOioSM predicts relatively large first and second generation scalar masses and 
smaller gaugino masses. An immediate consequence of such heavy first and second gen- 
eration sleptons is the suppression of the SUSY contribution to the anomalous magnetic 
moment of the muon. We find a^^^ < 3 x 10^^*^ (see Table). This is consistent with the 
most recent experimental and theoretical results at la if one uses r-based analysis pOj . 
However it is only consistent with an e"'"e~-based analysis at 2 a. 

Another important consequence of the MSOiqSM is the large value for tan (3 which 
leads to an enhanced branching ratio BR(i?s — > /u^ A*")- In addition, this is sensitive 
to the value of the CP odd Higgs mass tua Q, scaling as m~^. For = 500 GeV, 
the branching ratio satisfies 2 x 10^'' < BR(i?s /i+ /i^) < 5 x 10^'' for acceptable 
values of ^^^"^ ^"^^ ^ ^-^^ -^^ addition, we find that as mji increases, the region 
of parameter space consistent with WMAP data is forced to larger values of and 
smaller values of the light Higgs mass ruh- Hence, we find an upper bound on ttia ~ 1-3 
TeV consistent with the light Higgs mass bound mh > 114.4 GeV. For niA < 1-25 TeV, 
we find BR{Bs — > /i+^^) > 10^*^. Hence all acceptable regions of parameter space lead to 
observable rates for BR(i?s — > ^+ ^ )■ 

We update our predictions for the cross section for elastic neutralino-proton scattering 
due to scalar interactions function of the BR(i?s ^ /i"*" ) (Fig. P) and the 

neutralino mass (Fig. ^) for all regions satisfying the collider constraints, 0.094 < Cl^h"^ < 
0.129 and < 3. For comparison, in Fig. |7[ we also show the bounds from the present dark 

^^Where the upper limit is set by the recent DZero bound. 
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matter searches. Over the next two to five years the experimental sensitivity is expected 
to gradually improve by some three orders of magnitude. This will cover large parts of the 
predicted ranges of (7^^. Finally, we extend our previous analysis to include WIMP signals 
in indirect detection. 
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Table 1: The input parameters, ^ fits and Higgs and SUSY spectra for four representative points 
consistent with all available data. 
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